y=@(x)x^4-4*x^3-6*x^2-16*x+4;
c=[1,-4,-6,-16,4];
roots(c)
xmin=fminbnd(y,0,6)
polyval(c,xmin)
